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Quite general, analytical (both exact and approximate) forms for discrete 
probability distributions (PD's) that maximize Tsallis entropy for a fixed 
variance are here investigated. They apply, for instance, in a wide variety 
of scenarios in which the system is characterized by a series of discrete 
eigenstates of the Hamiltonian. Using these discrete PD's as "weights" 
leads to density operators of a rather general character. The present study 
allows one to vividly exhibit the effects of non-extensivity. Varying Tsallis' 
non-extensivity index q one is seen to pass from unstable to stable systems 
and even to unphysical situations of infinite energy. 

1 Introduction 

Tsallis' thermostatistics is today a new paradigm for statistics, with applica- 
tions to several scientific disciplines [21 EH EU E] • Notwithstanding its manifold 
applications, some details of the basic thermostatistical formalism remain un- 
explored. This is why analytical results are to be welcome, specially if they are, 
as the ones to be here investigated, of a very general nature. We will provide 
this type of results for discrete probability distributions of fixed variance that 
maximize Tsallis' entropy. 

Given a discrete probability distribution (DPD) p — {pk}, its Tsallis' informa- 
tion measure (or entropic form) is defined as 
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(1) 



It is a classical result that as q — > 1, Tsallis entropy reduces to Shannon entropy 

+00 

Hi(p) = - ^2 Pk^OgPk- 

fc— — 00 

Without loss of generality, we will here consider only centered random variables 
of fixed variance. 

The aim of this paper is to provide accurate estimates of i) the parameters 
of the DPD's and ii) their behavior. The maximizers of Tsallis' information 
measure under variance constraint in the continuous, multivariate case have 
been discussed in 0. 

Consider a quantum system whose eigenstates are characterized by a set of 
quantum numbers that we collectively denote with an integer, running index fc 
(Cf. Eq. QJ), that will of course also label the eigensolutions (\ipk), £fc) of the 
pertinent time-independent Schrcedinger equation 

W |Vfc> = £fc HAfe), (2) 

with H. the Hamiltonian. Let p^. = \ipk\ 2 be the probability of finding our system 
in the state |?/>fc). The mixed state 

p = X! Pk \i>k){i>k\, (3) 

fc 

commutes with the Hamiltonian by construction and represents thus a bona fide 
possible stationary state of the system. If we now find a physical quantity Z 
whose mean value is proportional to the variance, we can interpret p as the mixed 
state that maximizes Tsallis measure subject to the a-priori known expectation 
value of such a physical quantity. We will show below that, in these circum- 
stances, 



universal expressions can be given for the Pfc's, and thus for 



We discuss possible applications in the forthcoming section. Afterwards, after 
having introduced some definitions and notations, we characterize the discrete 
Tsallis maximizers for fixed variance in both the q < 1 and the q > 1 cases, 
and discuss thermal stability questions. We pass then to analyze extensions 
to the multivariate cases. For the sake of completeness, some of the proofs of 
assertions referenced to in what follows are given in Annex. 



2 Possible physical applications 

Several physical models can be adapted to the scenario described above (Cf. 
Eqs. (0 and J2J) The weights pk in (j^J) will be of the Tsallis-power law form. 
With these weights, p maximizes Tsallis' entropy subject to the constraint of a 
constant variance, which introduces a Lagrange multiplier that we will call f3. 
Given a physical quantity Z whose mean value is proportional to the variance, 
p is that mixed state which maximizes Tsallis measure subject to the a-priori 
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known (Z) — value. In the examples below Z is the system's energy E, but 
many other possibilities can be imagined. Thus, p will be the state maximizing 
Tsallis' H q for a fixed value of the expectation value U — (E) — Tr [H p] of 
the Hamiltonian. As a consequence, the multiplier (3 can be thought of as an 
"inverse temperature" T, that is, set (3 = l/(fcsT), with ks the Boltzmann 
constant. This is so because we are at liberty of imagining that U is kept 
constant because it is in contact with a heat reservoir 0. Of course, this is 
not necessarily the case, p exists by itself and is a legitimate stationary mixed 
state of our system. But we can think of (3 as either a "real" or an "equivalent" 
inverse temperature. 

Consider, for example, a system for which the energy spectrum consists of a 
denumerable set of N (N possibly infinite) energy levels labeled by a quantum 
number k with = Pk, so that all levels exhibit a degeneracy 

gk = 2 for k ^ 0; k > 0; and g = 1, (4) 

i.e., the sums in run from up to oo and each summand is multiplied by gk- 
Within the present framework, U 1 

U= (E) =J29kPkE k , (5) 

k=0 

becomes numerically equal to the DPD's variance, which by definition is fixed 
and assumedly known a priori. As just stated, we may think, if we wish, that 
our system is in contact with a heat reservoir, which fixes the mean energy, 
and that the associated Lagrange multiplier (3 can be assimilated to an inverse 
temperature T. Among many examples of such a scenario we mention here: 

• the planar rotor ; 8 j , where k is the magnetic quantum number correspond- 
ing to the azimuthal angle usually denoted by <p. The level- energies Ek 
are proportional to k 2 and 

Ek = Csk 2 ; Ce has dimension of energy. (6) 

We have Ce — h 2 /2Mj, with Mr the system's moment of inertia [§]. For 
simplicity's sake we take here Ce — 1, but retain, of course, its energy- 
units. 

• the three-dimensional rigid rotator [5], although k now means the quan- 
tum number L associated to orbital angular momentum (for large k, the 
spectrum of JBJ looks like that of the 3D rigid rotor) and, for all k > 0, 
9k = 1. 

1 Mean values linear in the PD are employed here. They are quite legitimate and were 
used by Tsallis in his seminal 1988 paper [2]. The Legendre structure of thermodynamics is 
definitely respected using them in conjunction with Tsallis' entropy 0. Recently, Ferri, Mar- 
tinez and Plastino have shown |4J that PD's obtained in this manner can be easily translated 
into PD's constructed via MaxEnt using q— expectation values evaluated a la Tsallis-Mendes- 
Plastino 5; , showing thereby that there is a one-to-one correspondence between the two types 
of PD. 
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• a vibrating string of length 21 with fixed ends whose energy eigenval- 
ues are Ek oc [kir/21] 2 . 

• the case of a particle (of mass m) in a box subject to periodic boundary 
conditions (at (— L, +L)) whose energy values are = [2n 2 H 2 k 2 / (mL 2 )] 

m 



3 Definitions and Notations 

In what follows, m is a positive real number and n a positive integer. Often, m 
will take the special form of an odd integer m = 2n + 1. The solutions to the 
problem of maximization of Tsallis entropy "under variance constraint" (here 
equivalent to "for fixed mean energy" ) will be called discrete Tsallis-probability 
laws (DTPL's) in this manuscript. Since (see preceding Section) fixed mean 
energy can be thought of as implying contact with a heat reservoir at the fixed 
temperature T = l/(ksP) we will call the pertinent Lagrange multipliers (3 
inverse temperatures in what follows. DTPL's are given by the following two 
theorems. 

Theorem 1 If h < q < 1, for ksT = 1/(3, with ks the Boltzmann constant, 
the discrete probability law p with zero mean and mean energy U ( equivalently, 
variance a 2 ) that maximizes Tsallis' entropy is defined as 2 

Pk = Pr{X = k} = f-\T)(l + Vfc e Z+, (7) 

where f q {T) is the partition function, given by 

is? p 

k=0 B 

and (3 = 1/fcgT is a real positive Lagrangian multiplier such that 



U = a 2 = y^gfcfc 2 Pfc- 



k=0 

In the limit T — > oo, our probability distribution converges to the classical dis- 
crete Gaussian distribution 

Notice the existence of a lower bound q = 1/3. Smaller q values are unphysical 
because for them the mean energy of the model (for the PD (7J|) becomes infinite. 
Additionally (see below) the system turns out to be thermodynamically unstable 
because the specific heat becomes negative. In this case where q < 1, we define 
the real positive parameter m that we need below by 

n.-l±£ (9, 



J the dependence of p on j3 = 1/fcsT is omitted for notational simplicity 
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Theorem 2 If q > 1, the discrete probability distribution p with zero mean and 
mean energy U ( equivalently, variance a 2 ) that maximizes Tsallis entropy is 
defined as 

Pk = Pr{x = k} = { tfPW-iSr)* Vfc e [0, Lfc B Tj - i] 

I otherwise 
where f q (T) is the partition function 

f q (T)= J2 ^-t-t^' 

fc=0 5 

[fceTj denotes the integer part of ksT whose inverse, namely, (3, is a real 
positive Lagrangian multiplier such that 



U = o- 2 = ^g k k 2 p k . 



k=0 

In the present case where q > 1, we define m as a real positive parameter related 
to q by 

»-!±i, (12) 

to be of use below. Note that this definition differs from definition © in the 
case q < 1. 

We do not give the proofs of theorems (1) and (2) here. That laws defined by 
JHJ and {IJl are the solutions of the Tsallis entropy maximization problem can 
be obtained by extending the results of the continuous case as presented in 
The cases i) q < 1 and ii) q > 1 differ essentially by the fact that the latter has 
a finite support. As far as we know, this is a novel way, in the q— literature, of 
dealing with the celebrated Tsallis-cut-off [5]. 

We remark on the fact that the forms J7J and (|l(Jfl coincide with the distribution- 
form given by Tsallis in his pioneer 1988 paper |2j (see also namely, 

p k = Z- 1 [l-(q-l)p*k 2 } 1 ^ 1 \ (13) 

that, for q — > 1, tends to the discrete Gaussian [2] 

Zf 1 cxp[-/Tfc 2 ]. (14) 

In the former case (Eq. (0) we have q < 1, while in the later (Eq. JIDJ) q > 1. 
Thus, in Eq. Q we have /? = (1- ^)/3*, while for 1(10)1 = (a- l)f3*. Of course, 
according to Q , fl9 > an d 03 , for q ^ 1 we have j3* — > /3. 
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4 Approximate treatments for q < 1 



4.1 The two regimes 

Let us recall that, for q < 1, m = (1 + <?)/(l — <z). Thus, m is large if q — ► 1 
(Boltzmann limit (BL)), while a small positive value of m entails q — > — 1. 
Careful inspection of the partition function f q (T) in the cases to = 3 <-> g = 1/2, 
to = 5 «-> g = 2/3, and to = 7 <-> 3 = 3/4 indicates that two regimes should be 
distinguished: 



the first regime, called "power law" regime, corresponds to the case T< 1 
(low temperatures) 



• the second regime called the "Student-t" regime, corresponds to the case 
T> 1 (high temperatures) 

Anticipating the results provided by the following theorems, we note that the 
cases "small T"/" large T" can be translated as joint range values of parameters 
U and to = (1 + q)/(l — q) : typically, the large T case corresponds to jointly 
large values of U and to, as expected, while the small case corresponds to jointly 
small values of these parameters. This is illustrated on the curve below (Fig. 1), 
where the large T property is translated as T > 100, and the small T property 
as T < 0.01. Note that the left bound n = 1/2 corresponds to to = 2n + 1 = 2, 
i.e., q — 1/3, and thus to an (unphysical) infinite U. For economy's sake we set 
herefrom 



a = v/A^T. (15) 
4.2 The power law regime 

In the power-law regime (T <C 1), a detailed characterization of the distribution 
can be obtained for all real values of to, as shown in the following theorem. 

Theorem 3 Assuming that T < 1 (a <C 1) and with to defined by pfy. the 
following approximations hold: 

1. for the partition function 

f q (T)^l + 2a m+1 am+l), 

2. for the probability law 



Pk 



-m — 1 



l +2a ™ + iC(m+l) n ' «/fc^0^ 
l+2a™ + 1 C(m+l) if k = 



and corresponds thus to a discrete power-law. It tends to a Gaussian for 



G 



3. for the mean energy 

C/~2a m+1 C(m-l). (17) 

4- If, moreover, m = 2n + 1 with n S N, n > 1, the characteristic function 
can be approximated as 

0m(M) - (l + 2a2«+2C(2n + 2)) (1 + 2 ° (2n + 2)! B2 ™+ 2 ^^' 

where B n (u) denotes the Bernoulli polynomial of order n. 
Proof, the proof is given in Annex 1. ■ 

The specific heat C = dU/dT is proportional to (m + 1). This is positive, and 
thus the system stable, for all ^ < q < 1. 

Notice that, according to l|16(l . only the ground state is populated at T = 0, as 
it should. 

The probability pk decreases exponentially with fc, being thus concentrated 
around its mean value EX = 0. As an example, if a = 0.1, the following table 
shows, for several values of q > i, the first values of pf.. 





k = 


k = ±1 


k = ±2 


k = ±3 




0.999 78 


9.997 8 x 10- 5 


6.248 6 x IO" 6 


1.234 3 x 10~ B 




0.999 98 


9.999 8 x 10~ b 


3.124 9 x 1Q-' 


4.1151 x 10~ s 




- 1.000 00 


1.000 00 x 10" B 


1.562 5 x 10~ 8 


1.3717x 10~ y 


1=7 


- 1.000 00 


1.000 00 x io-'' 


7.812 5 x 10~ 1U 


4.572 5 x IO" 11 



Tabic. 1: first values of pk for several values of q and a = 0.1 



4.3 The Student-t regime 

In the case of the Student-t regime (high temperatures), our main result writes 
as follows: 



Theorem 4 Ve > 0, 3oq > and 3uiq > such that if a > ao and 2 < m < mo 
then 

£5 fc 2 .™^ ^ r(f) 

2^fffc( 1 + -2) 2 - a V^7^I+T 

fc=0 i( - 2 

The proof of this result is given in Annex 2. The values of do and mo such that 
l|18|) holds are given in Fig. 2 for e = 10~ 9 . As a direct consequence of this 
result, we have the following explicit approximations, valid if, simultaneously, 
i) a > 1 (T > 1) and ii) q -> 1. 

2 

Proposition 5 I/a>l,m»l, and ^- ^> 1, men £/ie following approxima- 
tions hold with high accuracy: 



< e. 



(18) 
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for the partition function 

r(— ) 

f q (a) ~ a V^ r ^ m +iy 
/or i/ie probability law (remember that U = a 2 ) 

p ^ r( f) r(U^ (1 + ^^) r ^' (19) 



• lo/u'ie /or the mean energy 



U~^—. (20) 
m-2 K J 



Note that here 

(m-2)U = k B T = (21) 
a result that will be needed below. 

Proof, the expression of the partition function results from theorem fQ). Using 
a by-product of proof of theorem JjJ), the mean energy U reads writes 



C/ = a 2 (^-l)^a 2 (^=iIi-l) = 



/m JXf) m-2' 

r(=fi) 



^4s a consequence, a ~ yJ~U{m — 2) and i/ie Zaw writes 

T { m±l) e 
Pk T(f)T(i)^U(m-2y U(m-2)> 



Notice that in this case the specific heat is proportional to (m — 2) -1 . This is 
negative, and thus the system unstable, unless q > 1/3. In other words, the 
system is stable in the interval ]l/3, 1]. We remark also that, as a consequence 
of Ijl7|) and (|20|l . the mean energy U is an increasing function of a (and thus of 
T), as it should. 

Remark 6 A remarkable result here is that expression is exactly the sam- 
pled version - with the same partition function - of the continuous maximizer of 
the Tsallis entropy whose expression is recalled here 

f(x) = ; * ' (1 + TU Vx e R 

r(%)T(±)y/U(m-2y U{m-2)> 

This entails that, at very high T , one can approximate sums over discrete 
energy levels by an integral. In order to understand the situation, remember that 
the results holds for m 3> 1, which entails q ~ 1, i.e., the Boltzmann limit. In 
this limit, replacing sums by integrals is a commonplace text-book procedure. 



8 



5 The q > 1 instance 

5.1 A general result and its consequences 

Recall that in the case q > 1, m = (q + l)/(q — 1). The equivalent of theorem 
Q writes as follows: 

Theorem 7 Ve > 0, 3ao > and 3mo > such that if a > ao and m > mo 



J» y r(f + i) 

k=-(a-l) V 2 



< e. (22) 



Proof. The proof of this result is given in Annex 3. It essentially follows the 
steps of the proof of theorem Q). ■ 

We depict in Fig. 3 the area (north-east), delimited in the plane (o, m), that 
contains values of a and m for which approximation IL'l'l) holds, with e = 10 -10 . 

Corollary 8 For «> 1 (T>1) the mean energy U verifies 

a 2 

U ~ -2— (23) 
m + 2 v 7 

Proof. Denoting n — and f n (a) for f q (a), the normalization constant 

verifies the following recurrence 

+ a 1,2 + a 1,2 i + a 1.2 

fn+1 (a) = y: a - ^) n+1 =E( 1 4) n -iE fc2 (! - 

k——a k——a k——a 

= fn(a) - -n/n+i(a) 



that 



U = a 2 (l-l^). 

fn{a) 

Moreover, using the result of theorem TO, we /iat;e 
we obtain after some algebra 



U~-?—. (24) 
m + 2 v ' 



We see from the form of l|24|) that the specific heat C is proportional to l/(m+2). 
Thus, the is system stable for all g > 1. Also, since a 2 = ksT, we see that the 
mean energy U is an increasing function of T, as it should. 
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5.2 Convergence to a sampled Student-r law 

As a consequence of theorem |J7J| , the Tsallis maximizers can be approximated 
as follows, for a large enough (T large enough). 

Theorem 9 If a 1 then the following approximation holds 

Pk ~ ^ r( ^)V[7V^T2 ( ~ m^TW ■ (25) 

Proof. ^From we have 

a 2 ~ (m + 2)U. 
and the result follows from Theorem Q). ■ 
Note that here 

(m + 2)U = k B T = /3- 1 , (26) 

a result that will be needed below. We remark the similarity to the correspond- 
ing q < 1 probability law. 

Eq. (|25|l corresponds to the sampled version of a continuous Student-r distri- 
bution -which maximizes Tsallis entropy for fixed mean energy fTSJ - with the 
same partition function. 

5.3 Detailed results 

In the case where m is an odd integer, more detailed results can be obtained 
concerning the behavior of this probability law. 



Theorem 10 For a> 1, if m = 2n+ 1 is an odd integer with n > 2, then 

2 2 " +1 (n!) : 
(2n+l)! 



f m (a)= ,„ a + o(a- 2 ) (27) 



and Va > 1 , and if m = 3 



, n 4 a" 1 
h(a)=a--—. 



Moreover, the mean energy U verifies 



U=^-+o(a- 1 ). (28) 
m + 2 

Proof, the proof is given in Annex 4- ■ 

This quasi linearity of f m (a) in a for high temperatures is illustrated in Fig. 4, 
while the quadratic behavior of U vs. a is depicted in Fig. 5. 
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6 Convergence to the discrete normal distribu- 
tion 



We show now that, in the case where parameter m grows to infinity, both 
discrete Tsallis distributions corresponding to either q < 1 or q > 1 converge to 
the discrete normal distribution. 

Proposition 11 //«> 1 and if the mean energy U is fixed to a constant value, 
then, as m grows to infinity (i.e. q — ► 1 ), both discrete Tsallis probability dis- 
tributions Fw and $10)) converge to the maximizer of Shannon's discrete entropy, 
i.e., the discrete normal distribution. 

Proof, in the case q < 1 , the result follows by taking m — » +oo in fig)) with 

r(^±i) i 

hm — — = , 

m^+oo r(f)r(i) x /J7(m - 2) VWn 



and (Cf. JUPJ 



,. n k 2 .mil . k\ 

lu ? 1 + 777^ — ov 2 =ex P-^F7 
m^+oo (7 (m — 2) 2(7 



In the case q > 1, £/ie result follows similarly by taking m — ► +oo in i)£5j) urai/i 
lim 



r(f + i) _ i 



k 2 m-l fc 2 

lim (1 ; r)^~~ = exp( ). 

m^+oo V [7(to + 2) ; pv 2C/ ; 



This result shows that, for high enough temperatures, and if m — > oo (or 
equivalcntly g — > ±1) the distribution that maximizes Tsallis' entropy converges 
to the Boltzmann one, both for q < 1 and q > 1. Of course, for fixed T ^> 1 and 
(/ — > +oo, we obtain the uniform distribution, whereas for T<C 1 and q — > +oo 
we obtain the deterministic case pu = <5fc • 



7 Extension to the multivariate case 

In the last part of this communication we give a heuristic discussion of the 
multivariate discrete Tsallis laws with q < 1 . 

Theorem 12 if -A^ < Q < 1 fien f/ie probability Pk x ,...,k d = Pr{Xi = ki, ■■ ■ ,Xd — 
kd} defined by 
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k? 



=1 



Pk u ...,k d =/ m ,d(ai,---,a<i)(l + X)^) heZ+ 



maximizes Tsallis d-variate entropy 



H q ( P ) = ^j(l- £ Pl,... iki ) 

' k lt ...,k d 

for fixed mean energy 



ki 



ul 



The partition function f m ,d{o-i^ ■ ■ ■ , o-d) is defined as follows: 

d ^2 d 

f m ,d(ai,...,ad) = £ + 3_ 

ki,...,k d i=l 1 

In the case where Vi, a« 3> 1, £/ms function can be approximated as 
f m ,d{ai,...,a d ) ^ wm+d , 

1 V 2 / i=l 

Proof. The proof can be derived from that given for the continuous case in \@/. 
It is detailed in Annex 5. ■ 

The family of multivariate discrete Tsallis maximizers verify approximately, 
as in the continuous case, the property of invariance by marginalization, as 
described in the following theorem. 

Theorem 13 The d'—variate marginal law (d 1 < d— 1) writes 

k 2 ..+■>' 

Pk u ...,k d , = fm,d'(a-L,...,a d ,)(l +22^) 2 ■ (29) 

i=\ a i 

Proof, the proof is given in annex 5. ■ 

We remark here that property (|29|l is nothing but renormalizability: any subsys- 
tem of a system distributed according to a discrete Tsallis law is itself distributed 
according to a discrete Tsallis law. We note, too, that renormalization does not 
change the value of parameter m. 



8 Conclusions 

In this report, we have derived in analytic fashion some elementary properties of 
the discrete univariate and multivariate Tsallis laws (probability distributions 
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Pk running over a discrete index k) for both q < 1 and q > 1. Most of the 
ensuing results are given by approximations, but numerical simulations show 
that they can be regarded as very accurate ones. 

The discrete probability distributions pk define, for any Hamiltonian TL whose 
eigenvalue equation reads TL\ipk} — efclV'fc) mixed states of the form 



that commute with the Hamiltonian by construction and represent thus bona 
fide stationary states of the system defined by TL. These weights pk are of 
the Tsallis-power law form. With these weights, p maximizes Tsallis' entropy 
subject to variance constraint, which introduces a Lagrange multiplier (3. In 
this work we have considered situations for which the system's mean energy is 
proportional to the variance. 

Our present consideration allows one to nitidly appreciate the effects of non- 
extensivity, as, varying q, one passes from unstable to stable systems and even 
to unphysical situations of infinite energy. 
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FIGURE CAPTIONS _ 

Fig. 1: Plot of log 10 (VU) vs. [l/(q - 1)] for T = 100 (upper curve) and 
T = 0.01 (lower curve) that illustrates the regions of applicability of the Student - 
t and power law regimes. 

Fig. 2: Student-t regime for q < 1: the values of ao as a function of mo (see 
text). 

Fig. 3: The plane (a, m): values of a and m for which approximation \Tll 
holds, with e = 10~ 10 (see text). 

Fig. 4: q > 1 regime: quasi linearity of / m (a) in a for high temperatures. 
The four different curves correspond, respectively, to m=3,5, 7, 9, and 11, from 
top to bottom. Dots indicate the approximate value as given by Eq. (17). The 
continuous line represent exact values. 

Fig. 5: q > 1 regime: quadratic behavior of U vs. a. The three curves 
correspond, respectively, to m=3, 5, and 7, from top to bottom. Dots indicate 
the approximate value as given by Eq. (18). The continuous line represent exact 
values. 
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Annex: proofs 



For the sake of completeness we give here some of the proofs of assertions ref- 
erenced to in the text. 



8.1 Annex 1: proof of theorem (JIU) 

1. the normalizing factor writes 



— k\ m±1 k\ 



f q (a)= ^ (1 + ^)-^ =1 + 2^(1 

k — ~ OO A,'— 1 

+oo -,2 

- 1 + 2 J2 ( - — y Slt± = 1 + 2a m+1 C(™ + 1) 



fc=i 



2. the approximated expression of the distribution as given by l|16f) is a direct 
consequence of the preceding result 

3. We denote indifferently f m as f q . As a consequence of the preceding result 
and remarking that 



U = f~ 1 (a)a 2 J2 



+ oo i i k 2 +oo 

= /mW E v^ir-f^y E ttt^it 

= a 2 (%^-l) 
we obtain 

V / 2 «+i V l + 2a 2 "+ 2 C(2n + 2) y V 1 + 2a 2 ™+ 2 C(2n + 2) ' 

~ 2a 2n+2 (C(2n) - a 2 C(2n + 2))(1 - 2a 2n+2 C(2n + 2)) ~ 2a m+1 ((m - 1) 



4. if m = 2n + 1 with n£N, the characteristic function writes 

+°° 1,2 +°° u2 

4> x (u) = f m \a) J2 cos(ku)(l + -)-^ 1 =/-i( a )(l + 2^cos(^)(l + -)-^) 



* — ' ar 

k— — oo k—1 



fe=i 



But [E3 1.443] writes 



COs(fc^) (-!)»-! ( 2 ^) 2 " - 

fc=l 
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so that 



6 x(u )^/-i(a)(l + 2a-+ 1 ^^^; 

fc=i 

1 



1 + 2a m + 1 ((m + 1) 



(1 + 2a m+1 ( X) " ( 27r ) 2 " +2 R 2ra+2 ( — )) (0< U <2vr) 
V 2 (2n + 2)! 2n+2V 27r v ~ ~ 7 



8.2 Annex 2: proof of theorem flU) 

Denote the Gamma function as 

Using the Gamma integral, we have, with /x = 
+00 , 2 +00 „ +oc 



u 7 M (u)rfu. 



Consider an arbitrary value a 2 and divide the integral into two parts: 

- k2 v - k2 v f +OC - k2 v 

We will show that 

+ °° r +oo 



1 ..W 



can be made arbitrary small provided a is large enough, and for that choice of 

a, the absolute value of 7 2 — J^t^-oo Jo e~°? v 7n(v)dv — a V^~TJ]iy' can ^ e 
made arbitrary small by choosing fi small enough. 

In term I\, v > a 2 so that exp(— ^-w) < exp(— k 2 ) so that 



j-OO /- + 00 2 +OO /-+CSO M 

h= J2 e~^ v ltl {v)dv < J2 e_fe2 / 7/.(«)d« < 2 / 7/*(«)d« 

k=-oo Ja2 fc=-oo "' a2 "' a2 



since Ylk=-oa e fc2 — V 7 ^ — 1-772 6. As j a 2° ~fp{v)dv is the residue of the 
converging integral of a positive function, a proper choice of ao yields 

a> a < e/2. 
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The second term writes 



+ oc . a - fc2 

h = ^ / e~^ v ~f^(v)dv - 0^7? 
t._ „ Jo 



k— — oo 



I 2, e ** v Jft(v)dv ~ a V^ 



r(M-|) 
r(M) 

r(M-|) 



( / V e'^j^dv - / —^ ltl (v)dv) 
J k= -oo J o V« 



For the first difference, we use the fact that, for a > a\, 

lg,-fr_^i<. 

A;— — 00 



so that 

+ OC 



/" a J— IT 2 / — /" a 



This r.h.s. integral can be < \ as soon as a > ai. Choosing a > max(ai,a 2 ) 
yields 

/•a 2 +°c 2 <-a 2 /— 

1/ E e~^ v ^(v)dv- ^ lti {v)dv\<e/A. 

Jo u-~, Jo V v 



k— — OQ 

The second difference writes 



r 2 *v^ f x, ^ EQj-j) r- r 2 1 r w f +oc a ^ , w 

Jo V« r(;u) Jo Vv Jo Vv 

r(a--) r a2 

Now, for a > max(ao, ai) fixed, function /(/u) = a V^~TJJlf^ ln-^{ v )^ v ~ 1) 
is continuous in and such that 

Mm f(p) = 

so that 3/zo such that 

M<Mo =► |/(M)I <e/4. 
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8.3 Annex 3: proof of theorem 

Using the complex version of the Gamma integral as given in |131 3.382.6], we 
have 

p^= t m^ p r ,/ p \ _ dx 



2vr ' J_ x [fi + ixY 

if p > 0, Re(/3) > 0, Re(v) > 0. Choosing v - 1 = and /? = 1, we deduce 
that 



(1 - -^)~ = 2 e 1 ^ / ^rrefe 

a 2 2^ .Loo {1+ix) — 



Tf m+l \ r+ac i2 „l+ix 



2tt 7-oo (1+iat)' 



so that 



— fc=W r(^) r+°° ^ _ . 4 e i+- 



E (i-^ = iA ^/ E e M+1 cfa- 



fc=-(a-l) 30 fc=-(a-l) 

When a > 1 



E 



e -(i+«)^ 



so that 



, , s r(^) /-+ 00 a J^ e 1+ix , T(^)a^ f +ao e 1+ix 

f m (a) ~ — — / . m+1 dx 



2tt ./.^ v / T+«(1 + m) !; * 1 " 2tt + 

ii)aV^ 2vr a0Fr(2i±i) 



2tt r(f + i) r(f + i) ' 

8.4 Annex 4: proof of theorem fllOj) 

Denote n = and /„(a) = X)fc=-(a-i)(l — f^)" - ^ TO * s an integer then 
n is integer and f(n, a) can be expanded as follows 

a— 1 , 2 a — 1 n / \ ,2 n / \ a— 1 

/» = Ea - = E E (:) (-> = E (!) (-d p «- 2p E **• 

fc=0 fc=0p=0 p=0 k=0 

But it is well-known [T31 0.121] that, if p ^ 0, 

to 2 ^ +1 iV 1 j 2p+l 2 6 
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since B% = 0, and if p = 0, Ylk=o k 2p — a > so that 

n , «. a— 1 

= <• + £ Q <- i ) p »"' , (^T<' 2r+1 - r 2 ' + f + °<« 2p ")) 



(-I)p 0FI>+1) _ (2"n!) 2 



p— p— 1 x 7 p— 1 

Moreover, since 



E^^p+l 2 r(f + n) (2n+l)! 
by Euler's duplication formula, and as 

n>0 =► f^(-l)* = (l-l)»-l = -l 

and 



EC)(^--"5("; 1 )<- i ^-*-{" 1 oer 1 

we get finally, for n > 2, 

( 2 "™ ! ) 2 1 / -2s . / x 2(2"n!) 2 , _ 2 . 

f n {a) = H ho a 2 ) , f„(a) = a-± '— + o(a 2 ). 

3 y ' (2n + 1)! 2 y ' ' ^ v ' (2n + 1)! v ; 

and in the case n = 1, 

- . , 2 1a- 1 , , 4 a- 1 

/l(ffi)=a 3 + rT i/l(a) = B rT 



The partition function verifies the following recurrence 

w«o = £ a - ^) n+1 = E (i - -r -2 E fc2 (! - ^: 

' — ' a z z — ' a a z — ' n 

k— — a k— — a k— — a 



2/-, fn+lW- 



= fn(a) - -o/n(a) 

so that 

7n(o) 

Since Z7 > 0, this is an alternate proof that / ra (a) is decreasing in n. Moreover, 
using (|27l ). we obtain after some algebra 

C '=2^3 + <*■">■ 
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9 Annex 5: proof of theorems ( 1121) and ( 1131) 



We apply the same technique as in the univariate case, and represent the par- 
tition function as a mixture of discrete Gaussian using the Gamma integral as 
follows: 

d ^ 

f m ,d(ai,---,dd)= E ( 1 + E^2~) S_ 

ki,...,k d i=l 1 

= / e 1-1 '?7 m +d (y)dv = I e ~^j^ m+d (v)dv 

ki,...,kd ^ ^ ki,...,k d 

f+oc d _k'f f+oo d j— 

Jo i=i fc J o i=1 V« 



<i /- + oo ^ d 

7T2(TTai) / tl"?7 m-fj (v)dv = 7T^(TTl 

Jo 



Using the same approach, the d'-variate marginal law writes 



d , 2 



?*!,...,*„/= E P*i.-.*d = /m.d(«l--- a d) E ( 1 + E^) _ 













/mid( a l' ■ ' 


.,a<0 


E ca- 


i fr 2 

»=i * 


d 

E 




.,ad)(l 


4 — 1 ( 


m + d V— » 

2 E 


(1 + 

fed 



A Jfe?. 



' — ' a? 



with 



^ 2 = «?(i + E|) 



But, with to' = m + d' , 



e d+ e |r^-^( n ^ 

+ 1 ,...,fc d i=rf'+l 1 i=<i'+l 2 ^ 

d-d' "I r 1 l^-J ,„ \ ^ fcj , d-d' 

=* 2 ( n ^ w OTZ T a+E^) 2 

i=d' + l i( - 2 J i= l 
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so that 



Pk u ...,k A , 



d' 



(nti«,)r(f) 
r(*±£) 



i( - 2 i j=l 

2 



^(nti«0r(f 



■(1 



■(1 



_ h.2 
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